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Abstract 

The equations determining the energy density p of a gas of photons in ther- 
modynamic equilibrium with a spherical mass M at a non-zero temperature 
T s > is derived from Einstein's equations. It is found that for large r, 
p ~ 1/r 2 where the proportionality constant is a fundamental constant and 
is the same for all spherical masses at all temperatures. 
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§1. Introduction 

In the standard Schwarzchild solution [1] of Einstein's equations for 
a static, spherical geometry the spacetime outside of a sphere of mass M 
and radius R is taken to be a vacuum; empty and free of particles. This 
is certainly true if the sphere is at absolute zero temperature. When, how- 
ever, the sphere has a non-zero temperature, then we would expect, on a 
physically basis, a gas of thermal photons to be present. If the system is in 
thermodynamic equilibrium, then the spacetime surrounding the sphere will 
not be empty, but will instead be filled with blackbody radiation. As these 
thermal photons also have a certain non-zero energy, we would expect the 
presence of this blackbody radiation to also contribute towards determining 
the geometry of the spacetime. 

For any physically reasonable temperatures the energy density of the 
photons is very small in comparison to the mass density of the sphere and the 
presence of the photons is ignored in the usual Schwarzchild analysis. This 
approximation is certainly valid near the spherical body but what happens 
when one is very far away from the sphere? One should remember that 
the Schwarzchild solution is asymptotically flat. Due to the mass M of the 
body being confined within its radius R, as one goes further and further 
away from the sphere the affect of its mass M on the curvature of spacetime 
becomes less and less. The blackbody radiation, on the other hand, is 
unconfined. It extends from the surface of the sphere to fill the rest of the 
spacetime. Although we would expect the energy density of the photons 
to also decrease as one moves further and further away from the sphere, 
the important question is how fast it will do so. Let us, for the moment, 
consider a sphere of radius r > R and the total amount of energy within 
it. When r is near R, we would expect that the fraction €m of this energy 
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which is due to the mass M will be much greater than e 7 , the fraction of the 
total energy which is due to the photons. Consequently, we would expect 
the mass M to be the dominant factor in determining the geometry of the 
spacetime within r and would expect the Schwarzchild solutions to be valid 
in this region. As, however, r increases, €m decreases since the mass M of 
the body is fixed, while e 7 increases since the blackbody radiation extends 
throughout the spacetime. If the energy density of the photons decreases 
sufficiently rapidly so that e 7 <C 6m for all r, then the mass M will always be 
the dominant factor in determining the geometry of the spacetime. Although 
the Schwarzchild solution would be modified somewhat by the presence of 
the photons, we would not expect these modifications to be very drastic. In 
particular, the spacetime should still be asymptotically flat. If, on the other 
hand, the energy density does not decrease rapidly enough and e 7 ^> cm for 
r larger than some ro, then we would expect that in this region of spacetime 
it is the photons which will determine the spacetime geometry. In particular, 
we would expect the solutions to Einstein's equations in this regime to be 
very much different from the Schwarzchild solution. The spacetime may not 
even be asymptotically flat. 

In this paper we shall study some of the affects of non-zero temperatures 
on the static, spherically symmetric solutions of Einstein's equations. The 
system we shall be considering consists of a spherical body with a mass 
M and a radius R which, due to the sphere being at a temperature T s > 
0, is surrounded by blackbody radiation. The system as a whole will be 
assumed to be in thermodynamic equilibrium with the body serving as the 
heat reservior for the system. In particular, this means that the spherical 
body is assumed to be in thermodynamic equilibrium with the photons 
surrounding it. It is moreover assumed that the body has not undergone 

3 



complete gravitational collapse into a blackhole, is non-rotating, and is not 
electrically charged. Nor shall there be any other massive objects present in 
this spacetime. The only difference between this system and the one studied 
by Schwarzchild is the presence of a non-zero temperature for the sphere. 

Our aim in this paper is two fold. First, we shall derive a set of coupled 
differential equations which will determine the total energy density of the 
photons. As these are non-linear equations, we shall not be able to solve 
them analytically. We shall, nonetheless, be able to obtain both asymptotic 
r — > oo as well as r — > R solutions to them. Note, however, that we shall 
be determining the photon's total energy density in the gravitational field, 
and not its blackbody spectrum. Second, we shall show that the geometry 
of this spacetime differs drastically from the Schwarzchild geometry at large 
r. In fact with the presence of the photons the spacetime is no longer 
asymptotically flat. 

In order to avoid working with non-equilibrium systems, we have as- 
sumed that the mass M is in thermodynamic equilibrium with the blackbody 
radiation surrounding it. Unfortunately, even equilibrium quantum statisti- 
cal mechanics on curved spacetimes has yet to be satisfactorily formulated. 
The closest that we have come to a complete formalism is that given in [2] . 
For various reasons, however, it will be difficult, if not impossible, to analyze 
the system in the manner outlined therein and it is fortunate that all that we 
shall need is the Tolman-Oppenheimer-Volkoff equation for hydrostatic equi- 
librium. This, combined with the observation that the energy-momentum 
tensor of a gas of pure photons is traceless, shall be sufficient to determine 
the energy density almost uniquely. This method of deriving the energy 
density has the added advantage of not only taking into account the affects 
of the curvature of spacetime on the photons, but also the reciprocal affect 
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of the photons on the spacetime curvature. Photons are not treated as test 
particles in our analysis. In this, and other, ways our method differs from 
that discribed in [2]. The only difficulty that we shall encounter is when 
we try to identify the temperature of the system. Since we do not have an 
established formalism which will automatically do this for us, we shall, in 
the end, have to rely on other physical arguments to do so. 

In the formalism given in [2], as indeed in most treatments of ther- 
modynamics in general relativity, the temperature of the system is taken 
to be a constant throughout the spacetime. This is, it would seem to us, 
an oversimplification, for the following reason. For massless particles the 
temperature of the system at equilibrium may be interpreted physically as 
the most probable energy that any one particle in the statistical ensemble 
may have. (The case of massive particles is much more complicated and 
will not be considered here.) In a gravitational field this should include not 
only its kinetic energy, but also its gravitational energy as well. The two 
cannot be seperated covariantly. Indeed, it is known that the frequency of a 
photon, and thus its energy, when measured at different points in a gravita- 
tional field will either be "redshifted" or "blueshifted" with respect to one 
another. We would on this basis expect that temperature too should vary 
from point to point on the manifold. 

§2. The Hydrostatic Equation 

We begin with an iV dimensional manifold M with a metric which 
has a signiture of ( — ,+,+,+). Greek indices shall run from to N — 1 
and the summation convention is used throughout. It is further assumed 
that is static, meaning that there exists a timelike Killing vector £ M for 
the system. We shall also assume that the system contains only one heat 
reservior. 
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Next, let T Mi , be an energy momentum tensor operator defined on M. 
We shall, in a semi-classical approximation, treat g^ v as a background, clas- 
sical field. We next denote the thermodynamic average of T^ u by (T^). 
We shall not need a specific definition of this average, but rather that it 
satisfy a few basic properties that we would expect from any equilibrium 
thermodynamic average. First, it should be "time independent", meaning 
that 

^(V> = 0, (1) 

where £^ denotes the Lie derivative along the ^ direction. It is for this 
reason that we required M to have a timelike Killing vector. Physically, it 
means that the background field g^ v cannot change with respect to time the 
total energy contained in the matter fields so that the system as a whole 
can be in equilibrium. Second, we require that the average be anomaly-free 

V A (V) = (VaV), (2) 

where Va denotes the covariant derivative. Third, after the average is taken 
the energy momentum tensor should have the form (T^) — pu^u^ +p(g^ + 
UfiUu) were p and p are the proper energy density and pressure, respectively, 
and is a unit velocity vector which must lie in the £ M direction if (1) is 
to hold. We shall make this dependency explicit by writting 

{Tp V ) = -P^jf + P (g» v - ^f^J ■ (3) 

Then (1) requires that 

£ A V A p = , £ A V A p = 0, (4) 

meaning that p and p are functions of vectors lying in the N — 1 dimensional 
hypersurface perpendicular to ^ and/or £ 2 = Finally, since the 
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system is conserved, V M (T MZ ,) = and 

V M p+(p + p)^ = 0, (5) 

where |£| = a/~£ 2 and we have used (3) and Killing's equation 

+ V,£ M = 0, (6) 

in obtaining (5). This is the generalization of the Tolman-Oppenheimer- 
Volkoff [3] equation for hydrostatic equilibrium to general, static spacetimes 
and it reduces to the usual hydrostatic equation in the case of spherically 
symmetric spacetimes. 

We next consider a region of spacetime which contains only photons. 
In this region 

AttT^ = F mA F„ a - ^F Q/3 F^, (7) 

where is the field strength tensor. This energy-momentum tensor op- 
erator is traceless T£ = in four spacetime dimensions. Consequently, 
(Tff) = and for photons p = (N — l)p. The hydrostatic equation (5) is 
now trivial to solve yielding 

P = MP (8) 

where a is an arbitrary constant. The average energy momentum tensor for 
photons is thus given by 

(iv-i)ierV^ e 



= Tat T\T7Tn 9v ~ N J 7f ■ ( 9 ) 



To determine (T^„) completely one must first determine £ M for the manifold. 
As the presence of the photons will also affect the curvature of the spacetime, 
determining £^ ultimately involves solving Einstein's equations using (9) as 
the source term 

1 ° ( 9„-N%\ = ^R^ (10) 



(N-l) \£\ N V M J 8tt 
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where R^ v is the Ricci tensor and we are using geometrized units in which 
G = c = ks = h = 1. The term is absent since {T^J) is traceless. 

The task now is to interpret (8) physically. Let us consider, for the 
moment, the case of Minkowski spacetime and enclose the system in a very, 
very large box which is connected to a heat reservior at a fixed temperature. 
Killing's equation is now a simple partial differential equation and we may 
choose a coordinate system in which its solution for a timelike Killing vector 
is = (—ft , 0, 0, 0) where ft is a constant and the superscript / reminds 
us that this is the Minkowski spacetime. From (8) we find that in four 
dimensions, 

which, if we interpret 1/ ft as the temperature of the heat reservior, is just 
Boltzmann's law for photons. Then a = n 2 kg / (15h 3 c s ) is identified with 
the blackbody radiation constant. 

Using the Minkowski spacetime case as motivation, we shall tentatively 
identify the temperature T of the system as 

T = ip (12) 

in general, static spacetimes and see whether or not this will make physical 
sense. First, we note that although T does very with position, it is "time 
independent", namely 

£ t T = ^-£ M V M £ 2 = o, (13) 

as one would expect for a system in equilibrium. Second, we note that 
variations in T are due solely to the gravitational field. In fact, the temper- 
ature at various points of the manifold is related to one another by just the 
redshift factor, 

t (x) _(-e( X ')V\ (14) 



which is precisely what one would expect from time dilation and the fre- 
quency shift of photons in a gravitational field. Finally, we note that timelike 
Killing vectors £ in non-rotating systems are determined by the geometry 
of the spacetime only up to an overall constant, as can be seen explicitly in 
Killing's equation (6). Consequently, we have the freedom to attach an over- 
all constant to the Killing vector which we can then identify as the inverse 
temperature of the heat reservior. Indeed, from (14) we see that relative 
temperatures between two points on the manifold are determined solely by 
geometry and to determine an absolute temperature requires choosing a ref- 
erence point on the manifold from which we can measure all subsequent 
temperatures with respect to. The most natural reference point to choose 
is the heat reservior of the system and to measure all other temperatures 
based on its value. 

The procedure for determining the energy density of thermal photons in 
any static geometry is now clear. First we solve Killing's equations up to an 
overall constant for a timelike Killing vector ^ in terms of the componants 
<7 M z, of the metric of the spacetime. To determine the overall constant, we 
use as a boundary condition for £ M the temperature Th r of the heat reservior 
of the system by evaluating (12) at a point x^r on the surface of the heat 
reservior. It is required that the heat reservior have a constant temperature 
throughout its surface. Finally, Einstein's equations (10) are solved with 
the appropriate boundary conditions to determine g^ v . These boundary 
conditions are usually also given at the heat reservior. The energy density 
p of the photons and the geometry of the spacetime are thus determined. 

§3 Spherical Geometry 

We shall now attempt to solve (10) for a static, spherical geometry. 
Specifically, we shall consider a non-rotating spherical body of mass M, and 
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radius R which at its surface has a temperature T s . The spherical body will 
serve as the heat reservior for the system. The most general metric which 
is static and spherically symmetric has the form [3] 



ds z = -jdi £ + hdr z + r 2 d9 2 + r z sin 2 6d(j)\ 



(15) 



where / and h are functions of r only. Killing's equation is now straightfor- 
ward to solve giving £ M = (— Cfc/(r), 0, 0, 0), and £ 2 = —c\f{r) for a timelike 
Killing vector, is an arbitrary constant which is determined by evaluating 
(12) at the surface of the reservior: l/c^ = T s f{R) 1 / 2 . 

As for the boundary conditions for / and h, we choose them in the 
following manner. Take a point r > R and consider the amount of energy 
contained within a sphere of this radius which is due to the photons. Clearly, 
because of the presence of the photons, the geometry of the spacetime within 
r will be different from the Schwarzchild geometry. Let us, however, take 
r — > R + so that the amount of energy contained in the sphere due to the 
photons gradually decreases. Their affect on the geometry of spacetime 
must also do so correspondingly and just outside the body we would expect 
the geometry of the manifold to be the same as that of the Schwarzchild 
geometry. Consequently, we choose as our boundary conditions 



2M 

lim fir) = 1 



lim — - 

r^R+ h(r) 



1 - 



2M 



(16) 



Einstein's equations given by (10) are 
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where p = aT*f 2 (R) and the primes denote derivatives with respect to r. 
They may be reduced to two coupled, nonlinear differential equations 



327rpo 
P 



r + u 

rh V / h 



1 

rh 
h! 



, h 2 



1 / 1 



(18) 



It is doubtful that these equations can be solved analytically. There is, 
nonetheless, one general feature that we can determine from these equations 
alone. Since we require h > for all r > R, from (18) we find that / is a 
monotonically increasing function of r. Consequently, p = po/ f 2 is a mono- 
tonically decreasing function of r. The energy density of the photon gas is 
at its largest at the surface of the sphere and decreases monotonically as one 
goes further and further away from it. This is exactly what we would have 
expected physically. Notice also that since T = 1/|£| = T s [f(R)/f(r)] 1/2 , 
the temperature of the photon gas also decreases monotonically with r. 

We shall now obtain approximate solutions to (18) in the small and large 
r limits. We first consider the near field solutions when r is near R and write 
/«l-r, and h~ x w 1 — A for T, A^C 1. For the boundary conditions to 
be consistant with this approximation, we shall require 2M/R 1. Then 
ignoring terms quadratic in V and A, we find that 



h(r) 

f(r) » 1 



2M 


{■* 
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47rpo^ 3 
3M 



r 3 



r 
r 



1 - 



(19) 



Notice, however, that V and A increases quadratically with r and at some 
point r m they will no longer be valid. To estimate r m , we enforce the 
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condition that r(r) < 1 for r < r m . This gives 

[ rl ~ M, (20) 



47rpo 3 



as a determining equation for r m . The near field solutions (19) are valid as 
long as the total energy of the photons confined in a sphere of radius r is 
much less than the mass of the spherical body itself. 

As for the asymptotic, r — > oo solutions, we obtain them in the following 
manner. First we define 

'=5?' \ = '- 2K - < 21 » 

Then (18) may be written as 

dA 2A (2. ATjr ' 

A + 4K - 1 



1-2K V3 

^=A-K, (22) 

where y = log(r/ro) for some tq. These differential equations have a fix 
point at 

A a = K a = A (23) 

where the derivatives of A and K vanish. Then perturbing about this fix 
point, 

This matrix has eigenvalues 

so that the fix point (23) is stable. Physically, this means that no matter 
what initial conditions are chosen for A and K, both functions will eventu- 
ally flow to the fix point (23) at large enough r. We can see this explicitly 
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by solving (24) 



3 L . /r \ 3 / 4 . / V47 , r 

! 



47 / V47 , r 
cos I — log — 



1 . /V47 1 r 
+ i sml— log- 



(26) 



A, and ro are constants which require matching boundary conditions that are 
given at small r to determine. As the solutions to (22) for intermediate r are 
not known analytically, we are not able to do this explicitly. Nevertheless, 
numerical calculations, and a formal perturbative solution of (18) treating 
po as the perturbation indicates that ro ~ as l° n g as Vv^o > Trn ' 

We would therefore expect (26) to hold whenever r 3> l/y/p~o. Note also 
that solutions to (22) approach the fix point (23) very slowly; basicly as 

r -3/4_ 

In the very large r limit solutions to (22) asymptotically approaches 
f a = (567rpo/3) 1 / 2 r, and h a = 7/4 where we have used the subscript a to 

denote the asymptotic solutions. Thus the metric at large r is 

i 

ds 2 = - (^p) 2 r dt 2 + 7 -dr 2 + r 2 d6 2 + r 2 sin 2 ddtf . (27) 

and we can now see explicitly that this spacetime is not asymptotically flat. 
Next, we find that the asymptotic energy density is 

3c 4 

Pa = 567^' (28) 

where we have replaced the correct factors of c and G. At very large r, 
the energy density decreases as 1/r 2 with a proportionality constant which 
is an universal number and is independent of either the mass M or the 
temperature T s of the sphere. This is once again the consequence of (22) 
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being a non-linear differential equation and having a stable, non-zero fix 
point. The temperature of the photon gas in the asymptotic limit is then 



very large, one should remember that (29) is valid only when r is also quite 
large. As l p i ~ 10~ 33 cm, this ensures that ksT(r) will always be very much 
smaller than the Planck energy. In fact, (18) gaurentees that for r > R, 
T(r) < T s , the temperaure at the surface of the sphere. We should also 
mention that the asymptotic solutions are themselves solutions of (22) at 
any r, as can be seen explicitly. They do not, unfortunately, satisfy the 
correct boundary conditions. 

We can also calculate the total average energy of the photons 



from (18). Since m ~ r for large r, we do not expect this E to be finite. 
It should instead diverge linearly with r. This divergence is much milder, 
however, than in the case of flat spacetime where the total average energy 
diverges as the volume of the system. 
t§4. Discussion 

The spacetime outside of a sphere with temperature T s can thus be 
divided into three regions. In the near field region, r < r m and the solutions 
to Einstein's equations are given by (19). The geometry of the spacetime 
in this region is dominated by the mass M at r = and the presence of 
the photons will not have an significant affect on it. In the intermediate 
field region, r m < r < 1/ p Q and the total energy contained in the thermal 




(29) 



where m p i = (hc/G) 1 / 2 is the Planck mass and l p i = m p iG/c 2 is the Planck 
length and we have explicitly used a = n 2 k% / (lhTC 1 c 3 ) . Although m p i is 




(30) 
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photons is now comparable to the mass of the sphere. Both the photons 
and the mass M together will determine the geometry of the spacetime. In 
the far field region, r l/y/p Q and the asymptotic solutions (26) are now 
valid. It is now the photons which are dominant over the mass M. 

We have in this paper considered only systems which are in thermo- 
dynamic equilirbium. In particular, this means that the sphere must be in 
thermodynamic equilibrium with the gas of photons surrounding it. As the 
spacetime that we have been considering consisted of only the mass M and 
the photons, this is equivilant to saying the sphere must be in thermody- 
namic equilibrium with the rest of its universe. Since our universe is filled 
with with the cosmic microwave background radiation which is at a temper- 
ature of ~ 3° K, for a physical body to be in thermodynamic equilibrium 
with the rest of our universe it must also be at a compareable temperature. 
There are very few actual physical bodies which are at such low tempera- 
tures. Moreover, because the mass M is in equilibrium with the photons 
surrounding it, the amount of energy radiating away from the sphere must 
exactly be balanced by the amount of energy impacting on the sphere by 
the photon gas surrounding it. It is for this reason that the energy density is 
dependent only on the geometry of the spacetime, and is why the definition 
of the temperature (12) makes sense. It is also the reason why the intensity 
of the emitted radiation from the sphere does not have the characteristic 
1/r 2 behavior as one would naively expect. 
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